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As a common experience from basic courses in Classical Mechanics, for some mechanical 

^ . systems the equations of motion can be solved up to quadratures, while others persist to 

! deny such access. This experience can be formalized and leads to the notion of an inte- 

•/^ \ grable system. For a Hamiltonian system with n degrees of freedom, one requires to have 

^ I at least n functions on phase space, Hj, j = 1, . . . ,n, which are in involution meaning 

^ ■ that the Poisson brackets {Hi,Hj} = for i, j = 1, . . . ,n, see |T] for details. Hi, say, is 

O ' the system's Hamiltonian. Then the manifold {Hj = cj, j = 1, . . . ,n} has the structure 

^ . of an n-torus and the motion is characterized by at most n frequencies. Hence, up to 

^ i deformation, the motion looks like the well-known Lissajous figures. 

I The text book example is the motion of a particle subject to a central potential. More 

rS ■ spectacular is the observation that integrability persists for particular systems with a 

■ large number of degrees of freedom, which first surfaced indirectly through the discovery 



of solitary wave solutions by N.J. Zabusky and M.D. Kruskal [2J for the Korteweg-de-Vries 
equation in one dimension and for a chain of nonlinear coupled oscillators by M. Toda 
[3]. A very rich field ensued. In the following my focus will be on the aspect of many 
interacting components. 

Naturally one may wonder how integrability survives under quantization. An Hamilto- 
nian operator, H, allows for many commuting operators. Thus a simple minded extension 
from the classical case will not do and there seems to be no generally agreed upon defi- 
nition of quantum integrability. On the other side there are clear signatures to identify a 
quantum integrable system (once it is found), to name only a few: Bethe ansatz, Yang- 
Baxter equation, and factorized S'-matrix. 

From the perspective of statistical mechanics it is also a natural issue to understand 
whether and how integrability extends to stochastic systems. To have one distinction 
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very clear, many systems of 2D equilibrium statistical mechanics are integrable, the cor- 
respondence being related to the fact that the transfer matrix has a structure akin to a 
quantum integrable system. In contrast, here I discuss stochastic time evolutions mod- 
eled as a Markov process, either diffusion or jump. As a linear operator, the generator, 
L, of the Markov process has possibly some structural similarity to —H, hence it seems 
reasonable to expect a corresponding version of integrability. On the other side, e^* is 
already the normalized transition probability; there are no probability amplitudes, the 
partition function equals 1, and the largest real part of the eigenvalues is 0. 

With R. Dobrushin the Russian probability school pioneered the many component 
aspect. Integrability is usually first associated with R. Glauber's exact solution of the 
one- dimensional stochastic Ising model [1]. This solution is based on what is now called 
duality, a concept introduced and generalized to other systems by F. Spitzer in the very 
influential article [5]. The dual description is here in terms of evolution equations for 
the time-dependent correlations functions, which decouple for integrable systems. An 
example is the symmetric simple exclusion process on the one- dimensional lattice Z. In 
this model there is at most one particle per site and, under this restriction, particles 
perform independently nearest neighbor symmetric random walks. The generator L equals 
—H with H the Hamiltonian of the ferromagnetic Heisenberg chain. (In this case, duality 
holds in arbitrary dimension and also for longer ranged symmetric jumps.) 

On the level of duality, none of the signatures known for quantum integrability make 
their appearance. This situation changes drastically as we turn to the asymmetric version 
of the simple exclusion process, ASEP (now ID and n.n. do matter). A particle at site 
i jumps to site j + 1 with rate p and to site j — 1 with rate g, g -|- p = 1, provided the 
destination site happens to be empty. The symmetric case corresponds to q = p = |. The 
generator can be written in the notation of quantum spin chains. If cr| = 1 means site j 
is occupied by a particle, then 

^ = i E (^^- ■ ^^+1 - 1 + 2^(p - - ^M+i)) • (1) 

Note that L is not symmetric. All eigenvalues are in the open left hand plane except for 
0. On a ring with a fixed number of particles, the unique invariant measure is the uniform 
distribution. The other eigenvectors are determined through the Bethe ansatz [6]. Much 
more powerful is the Bethe ansatz inspired expression for the transition probability e"^* 
discovered by C. Tracy an H. Widom [7]. Their expression is still extremely complicated 
and to simplify further one has to specify some initial conditions. A widely studied choice 
is the initial step, for which the half lattice {j < 0} is empty and {j > 1} occupied. For 
q > p Tracy and Widom write a Fredholm determinant for the probability distribution of 
Xj{t), the position of the j-th particle at time t. Much earlier K. Johansson found a 
distinct Fredholm determinant for a related quantity in the totally asymmetric limit q = 1 
(TASEP). Both results serve as the stepping stone for an intricate asymptotic analysis 
eventually arriving at objects familiar from random matrix theory. 

Very recently one accomplished to cross the border from discrete jump processes to a 
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particular stochastic PDE, which reads 

and is the ID version of the equation first proposed by Kardar, Parisi, and Zhang as 
a model for growing interfaces. Here h{x,t) is viewed as a height function and W{x,t) is 
white noise in space-time. Integrability is seen most explicitly for the sharp wedge initial 
condition, 

1 

h{x,0) = --\x\ with 5^0, (3) 



which, at 6 = 1, should be understood as the analogue of the once integrated initial step. 
(|2]) together with (|3]) looks very singular, and it is. For smooth initial data the solution 
is constructed by L. Bertini and G. Giacomin fTO] and for the sharp wedge in 
The KPZ equation turns linear through the Cole-Hopf transformation 

Z = e\ (4) 

Then 

|z = l^Z + iyZ, Z{x,0) = 6{x), (5) 

from which one concludes that the exponential moments of h are linked to the attractive 
5-Bose gas in one dimension, which is a quantum integrable system solvable through the 
Bethe ansatz. For example, for ([2]) together with ([3]) it holds 

E(Z(0,t)") = (0|e"*^"|0) (6) 

with Hn the n particle attractive Lieb-Liniger hamiltonian, 

1 " (9^ 1 " 

j=l j i¥'j=^ 

and |0) the state where all n quantum particles are at 0. Unfortunately, the moments 
in ([6]) increase as exp(r;,^), which makes a rigorous control difficult. But replica schemes 
have been employed and yield fascinating results [121 [ISl 111 [13 HB] . 

Currently the integrability of the KPZ equation can be deduced only indirectly by 
taking a continuum limit of the asymmetric simple exclusion process, where the lattice 
spacing is e, the time scale £~^, and the asymmetry q — p = y/e with e ^ 1. To give an 
impression, I record the generating function for the height at the origin at time t, 

E(exp [ - e-^e'^W+^^Z^^)]) = det(l - Poi^s,tPo) • (8) 

Here the determinant is in L^(M), Pq projects onto [0, oo), and Ks^t is an operator with 
integral kernel 

Ks,t{x,y) = [ (l + e-(*/2)'''"+^)"'Ai(a; + A)Ai(y + A)dA (9) 
Jr 
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with Ai the Airy function. PoKg^tPo is of trace-class. For large t, h(t) = — 1/24 + 
where the random amplitude ^ is Tracy- Widom distributed, just as is the largest 
eigenvalue of a GUE random matrix in the large N limit. (j8]) together with ([9]) was 
obtained independently in [HI [TTl [18]. In this context the introductory review [19j is 
highly recommended with some complementary information provided in pOj. 

The integrability of the KPZ equation triggered further advances. One interesting 
direction is to consider discretized versions of the stochastic heat equation ([5]). Somewhat 
unexpectedly, the completely asymmetric discretization turns out to be more tractable 
and one starts from the equations of motion 

j^Z.it) = - z,it) + (|6,(t))Z,(t) , Z,(0) = Sjo , J e Z , (10) 

where {bj{t),j G Z} is a collection of independent, standard Brownian motions. N. 
O'Connell [21] established a close connection between logZ„(t) and the last particle in 
the open quantum Toda chain of n sites. Very recently A. Borodin and I. Corwin [22] 
explain how Macdonald functions enter the picture. They are the eigenfunctions of the 
commuting set of Macdonald operators. In the future, for sure, the interface between 
stochastic and quantum integrability will be further elucidated. 

While we emphasized the notion of integrability, let me add as a fairly extended foot- 
note that the predictions based on the exact solutions have been confirmed recently in 
spectacular experiments [23j, see also the expository article [24J. Of course, physical sys- 
tems are much more complex than simple models as the TASEP. But on a large space-time 
scale microscopic details hardly matter, except for generic properties, like the condition 
of a sufficiently local update rule. In fact, such universal behavior can be proved for 
the integrable models discussed, but it should hold at much greater generality, including 
physical systems. In the experiment p3] one studies droplet growth in a thin film of 
turbulent liquid crystal. The film thickness is 12 /zm, while the droplet grows laterally to 
a size of several mm. The droplet consists of the stable DSM2 phase and is embedded in 
the metastable DSMl phase. Hence the interface is a line and it advances through nu- 
cleation events where the stable phase is created out of the metastable one. On average, 
the solution to the KPZ equation with sharp wedge initial data has a parabolic profile 
which self-similarly widens linearly in t and thus models one section of the droplet. By the 
physical conditions, the droplet growth is isotropic guaranteeing that the non-universal 
coefficients do not depend on the direction of growth, which is the basis for high precision 
sampling of entire probability density functions. In fact, the GUE Tracy- Widom distri- 
bution for the height fiuctuations is confirmed with accuracy. It is also observed that 
for fiat initial conditions, h{x, 0) = 0, the height fiuctuations switch from GUE to GOE 
statistics, implying that some features of the initial conditions are still visible in the large 
scale universal limit. 
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